A PREDICTION OF OBSERVABLE ROTATION 
IN THE ICM OF ABELL 3266 



KURT ROETTIGER 

Department of Physics and Astronomy 
University of Missouri- Columbia 
Columbia, MO 65211 
E-mail: kroett@hades.physics.missouri.edu 



RICARDO FLORES 

Dept. of Physics and Astronomy 
University of Missouri-St. Louis 

St. Louis, MO 63121-4499 
E-mail: Ricardo .Flores@umsl. edu 



Accepted for publication in the Astrophysical Journal 



ABSTRACT 

We present a numerical Hydro+N-body model of A3266 whose X-ray surface bright- 
ness, temperature distribution, and galaxy spatial and velocity distribution data are 
consistent with the A3266 data. The model is an old (~3 Gyr), off-axis merger having a 
mass ratio of ~2.5:1. The less massive subcluster in the model is moving on a trajectory 
from southwest to northeast passing on the western side of the dominant cluster while 
moving into the plane of the sky at ~45 degrees. Off-axis mergers such as this one are 
an effective mechanism for transferring angular momentum to the intracluster medium 
(ICM), making possible a large scale rotation of the ICM. We demonstrate here that the 
ICM rotation predicted by our fully 3-dimensional model of A3266 is observable with 
current technology. As an example, we present simulated observations assuming the 
capabilities of the high resolution X-ray spectrometer (XRS) which was to have flown 
on Astro-E. 
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1. Introduction 

A3266 is a nearby (z=0.059; Quintana, Ramirez 
& Way 1996), X-ray luminous cluster which exhibits 
optical and X-ray substructure. Two models have re- 
cently been proposed to explain the substructure in 
this cluster. Flores, Quintana, & Way (1999) have 
proposed that A3266 experienced a major merger 
into the plane of the sky while Henriksen, Donnelly, 
& Davis (1999) have proposed a minor merger par- 
allel to the plane of the sky. Using a large sam- 
ple of galaxy redshifts (387 galaxies), Quintana et 
al. (1996) suggested that A3266 might have experi- 
enced a merger 1-2 Gyr ago. Flores et al. (1999) 
found support for this interpretation using a simple 
N-body model. They noted an enhancement of galax- 
ies north of the X-ray core similar, both visually and 
statistically, to the N-body particle spray found in 
their numerical simulations. They also noted a simi- 
lar enhancement of emission-line galaxies in the same 
region. It has long been suggested that galaxies pass- 
ing through cluster cores could be spectroscopically 
altered (Dressier & Gunn 1983), although it now ap- 
pears that this would be mostly due to the tidal force 
rather than ram pressure by the ICM (e.g., Moore et 
al. 1996; Bekki 1999; Fujita et al. 1999). Burns et 
al. (1994) have suggested that the Coma cluster E+A 
galaxies distributed in the core and SW toward the 
NGC 4839 group are the result of a burst of star for- 
mation induced by a merger about 2 Gyr ago, which 
appears consistent with their starburst age (Caldwell 
et al. 1996). Similarly, the emission- line galaxies in 
A3266 could be the tail end of the disrupted less mas- 
sive cluster in this model. 

Evidence of recent dynamical evolution is also ap- 
parent in the X-ray surface brightness (XSB; Fig. |]) 
and temperature distributions. For example, the XSB 
was shown to exhibit a systematic centroid shift by 
Mohr, Fabricant & Geller (1993). Also, the XSB ex- 
hibits changing ellipticity and isophotal twisting be- 
tween 4' and 8' (Mohr et al. 1993) as well as a 
large (^500 kpc) core radius (Mohr, Mathiesen & 
Evrard 1999), as would be expected in the case of 
a recent merger (Roettiger et al. 1996). Peres et 
al. (1998) find no evidence of a cooling flow. Several 
researchers (e.g., McGlynn & Fabian 1984) have sug- 
gested that mergers will disrupt cooling flows. Gomez 
et al. (1999) have demonstrated using numerical sim- 
ulations that the time scale for re-establishment of 
the disrupted cooling flow in the post-merger envi- 



ronment can be greater than several billion years, de- 
pending on the initial cooling flow and merger param- 
eters. Markevitch et al. (1998; hereafter MFSV98) 
and Henriksen et al. (1999; hereafter HDD99) have 
produced temperature maps based on ASCA data 
that show significant temperature variations across 
the cluster. MFSV98 shows a radially decreasing tem- 
perature profile ranging from 12 keV in the central 
3' to ~6 keV at radii greater than 8' . De Grandi 
& Molendi (1999) find a similar radial temperature 
gradient using BeppoSAX data. The HDD99 temper- 
ature map exhibits a comparable range in ICM tem- 
peratures with a gradient increasing from NE to SW 
across the cluster. 

In this paper, we extend the N-body model of Flo- 
res et al. (1999) by including the hydrodynamics 
of the ICM. We then demonstrate using a fully 3- 
dimensional numerical model that the current A3266 
data are consistent with an old off-axis merger oc- 
curring largely into the plane of the sky. Off-axis 
mergers are a natural consequence of large-scale tidal 
torques, the latter being a generic feature of hierar- 
chical clustering (Peebles 1969). The model of Flores 
et al. (1999) resulted in an off-axis merger as a re- 
sult of the global angular momentum imposed at the 
protocluster stage, which in terms of the standard 
dimensionless angular momentum A (Peebles 1969) 
corresponded to A = 0.07. This amount is consistent 
with tidal torquing and is expected to be largely inde- 
pendent of mass (Barnes & Efstathiou 1987). There- 
fore the characteristics of the merger are not sensi- 
tively dependent on extending the region that was 
simulated around the cluster. In sufficient quantity, 
angular momentum can significantly alter the inter- 
nal structure of clusters, which can then influence our 
interpretation of other cluster observations. As an 
example, numerical simulations (Inagaki et al. 1995; 
Roettiger et al. 1997) have shown that the shape of 
clusters (e.g., oblateness being a consequence of ro- 
tation) can have significant systematic effects on de- 
terminations of H based on the Sunyacv-Zcldovich 
effect (see Birkinshaw (1999) for a review). 

We describe our model and make direct compar- 
isons to the data in §§. In §||, we present detailed 
models of proposed Astro- E observations based on 
the line-of-sight (LOS) ICM density, temperature and 
velocity structure provided by the numerical model. 
Section ^ is a summary of our results. We assume 
H =70 km s _1 Mpc -1 when scaling the simulation 
to physical dimensions. 
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2. A Numerical Model of A3266. 

We have created a numerical model of A3266 using 
the same technique that we have employed in several 
previous models of specific Abell clusters (e.g., A2256 
Roettiger et al. 1995; A754, Roettiger et al. 1998; 
A3667, Roettiger, Burns & Stone 1999). Within the 
framework of idealized initial conditions, we survey 
merger parameter space (mass ratios, impact param- 
eters, gas content, etc.) using a fully 3-dimensional 
Hydro/N-body code based on the Piecewise-Parabolic 
Method (PPM; Colella & Woodward 1984) and a 
Particle-Mesh (PM) N-body code. We then attempt 
to maximize agreement between synthetic observa- 
tions of the simulation and various cluster observ- 
ables (X-ray surface brightness, X-ray temperature 
distribution, galaxy spatial and velocity distributions 
etc.) in order to constrain not only the merger pa- 
rameters, but also the epoch and viewing geome- 
try of the merger. Here, we have modeled A3266 
as an off-axis merger between a primary cluster of 
~T.lxl0 15 M© and a secondary of ~5xl0 14 M Q in 
which closest approach occurred approximately 3 Gyr 
ago. The secondary cluster, moving southwest to 
northeast, passed to the west of the primary clus- 
ter's core at a distance of ^230 kpc with a velocity of 
^2500 km s~ x . The trajectory is believed to be into 
the plane of the sky at an angle of ~45°. 

Figure || shows the synthetic XSB image generated 
from the model. Like the ROSAT PSPC image (Fig. 
Q), the simulated image shows a generally spherical 
distribution at large radii with significant isophotal 
twisting near the X-ray core which is elongated NE 
to SW. The orientation of the X-ray cores (both simu- 
lated and observed) are not well-aligned with the pro- 
jected mass distribution derived from gravitational 
lensing (Joffre et al. 1999) further indicating that 
the cluster is not fully relaxed. It should be noted 
that the resolution of the numerical simulations (~75 
kpc or ~4 zones core radius) is significantly less than 
the resolution of the ROSAT image (15" /pixel or 
~20/if x kpc). 

Also included in Fig. || is a sampling of ~300 N- 
body particles (<1% of total) from both the primary 
(o) and secondary (+) clusters. The excess of sec- 
ondary cluster particles to the north of the X-ray 
core accounts for the galaxy excess as well as the 
distribution of emission-line galaxies noted by Flo- 
res et al. (1999). The observed galaxy velocity dis- 
tribution is indistinguishable from Gaussian within 



the central region. The observed skewness and kur- 
tosis within the central 1° x 1° field are 0.024 and 
0.14, respectively, while the N-body particle values 
are 0.051±0.15 and 0.13±0.26. On a larger scale 
(-2°), Flores et al. (1999) find the observed skew- 
ness and kurtosis to be 0.106 and 0.341, respectively. 
The global N-body velocity dispersion is 905±40 km 
S _1 , while Quintana ct al. (1996) observe a global 
galaxy velocity dispersion of 1085±51 km s _1 . The 
discrepancy here is possibly due to choosing too low 
a value for the initial cluster /3-parametcr. 

A3266 contains a central dumb-bell galaxy with 
a velocity separation of ^400 km s -1 . Quintana et 
al. (1996) suggest that this is consistent with the 
merger geometry proposed here. Although consistent 
with a younger merger that is nearly in the plane of 
the sky, it is also consistent with an older merger at 
any projection since the time scales for dynamical fric- 
tion to act on the dominant galaxy are long (Kravtsov 
& Klypin 1998). The velocity separation noted here 
should be contrasted with the 2635km s _1 separation 
observed between dominant galaxies in the proposed 
young merger A2255 (Burns et al. 1995). 

It has been suggested that radio source morphol- 
ogy may give clues to the ICM dynamics. The ra- 
dio emitting plasma is believed to be 10 to 100 times 
less dense than the surrounding thermal gas making 
it susceptible to pressure gradients and flows within 
the ICM (e.g., Burns 1998). HDD99 used the mor- 
phology of two radio sources to support their model, 
and depending on the exact location of these sources 
within the cluster, they are also consistent with the 
model presented here. We should comment however 
that considerable caution must be used when inter- 
preting radio source morphology in this context. Both 
sources in question are located SW of the cluster core 
(see Jones & McAdam 1992). One appears to be a 
head-tail source while the other is identified as a pos- 
sible Wide- Angle Tailed radio source (WAT). WATs 
have been used as indicators of bulk flows in clusters, 
because they have traditionally been associated with 
central dominant galaxies which are presumed to be 
at rest in the cluster's gravitational potential min- 
imum making knowable its exact location (and rel- 
ative velocity) within the cluster. Consequently, any 
bending of the WAT tails is the result of ICM dynam- 
ics rather than motion of the host galaxy. However, 
both of these radio sources are associated with galax- 
ies having velocities significantly different (AV^800 
km s _1 ) from the mean cluster velocity indicating 



4 



that they may have considerable velocities of their 
own and may even be foreground or background ob- 
jects. 

Figure [3] shows the projected, emission-weighted 
temperature map overlaid with the XSB contours. 
Often the X-ray temperature distribution provides 
the strongest constraints on the merger parameters. 
Two X-ray temperature maps have been published 
recently based on ASCA data (MFSV98; HDD99). 
Although largely consistent, they do differ system- 
atically. Both maps show a similar degree of tem- 
perature variation (^6-12 kev) within the cluster. 
However, the MFSV98 map shows a hot core with 
a radially decreasing temperature profile while the 
HDD99 map shows more of a temperature gradient 
across the cluster. The core is not the hottest region 
in the HDD99 map. For this reason, we have per- 
formed a region-by-region comparison of our model 
with both published temperature maps. Figure ^a 
is a comparison between our model and the corre- 
sponding regions in the MFSV98 map. Figure |Jb 
compares the model temperatures within regions de- 
fined by HDD99 (see Fig. |for region definition) . Our 
model agrees within the 90% confidence intervals of 
all but one region defined by MFSV98. In an absolute 
sense, the agreement with the HDD99 map is just as 
good although their quoted uncertainties are signifi- 
cantly smaller than those quoted by MFSV98. The 
most significant discrepancy between our model and 
HDD99 is in region 5 which includes the cluster core. 
Although our model agrees quite well with MFSV98 
in this region (Region 1; Fig. ^a), the HDD99 map 
indicates a cooler core. If correct, this could be an 
indication that cooling may have started to influence 
the core (where cooling times are the shortest) be- 
cause the merger is relatively old. Our current model 
does not include the effects of radiative cooling. 

A3266 has recently been observed in the hard X- 
ray band (15-50 keV) by BeppoSAX (De Grandi & 
Molendi 1999). They found no evidence of a hard X- 
ray component nor is there currently any evidence of 
diffuse radio emission. This result is consistent with 
an old merger interpretation since shocks present in 
young mergers might be expected to accelerate rela- 
tivistic particles which may produce observable dif- 
fuse radio emission via synchrotron and hard X-rays 
via inverse Compton scattering (see Sarazin 1999 for 
a recent review). 

Given the model presented here, what can we say 
about the gas dynamics in A3266? First, A3266 



appears to be an old off-axis merger in which the 
angular momentum of the initial clusters is due to 
tidal torquing (Flores et al. 1999). Previous off-axis 
merger simulations (Roettiger et al. 1998; Ricker 
1998) have shown that mergers can transfer signifi- 
cant angular momentum to the ICM and that this 
angular momentum can be long-lived. In fact, be- 
cause of the time and distance scales involved, it takes 
a considerable period (several billion years) for full 
rotation to develop. To understand this, consider a 
mean rotational velocity, v ro t=1000 km s _1 , at a ra- 
dius, r =500 kpc. The time required to complete cir- 
culation about the core is then simply t — lixr /v ro t or 
~3 Gyr. In this model, full rotation occurs between 
2.5 and 3 Gyr after core passage and persists beyond 
4 Gyr. Of course this is in the absence of a second sig- 
nificant merger event which could potentially disrupt 
the circulation. 

Figure || shows the gas velocities within an east- 
west plane taken along the observer's LOS and through 
the cluster core. Note that even at this late stage of 
the merger, there are still significant gas velocities 
(>500 km s _1 ) and that full rotation has been es- 
tablished. Further note, that shocks generated early 
in the merger have now dissipated while extremes in 
the temperature distribution have not. As we have 
pointed out previously (Roettiger et al. 1996), sub- 
structure of this type can persist well beyond the 
canonical sound crossing time. Thus, it may signif- 
icantly influence estimates of cosmological parame- 
ters based on the frequency of substructure in clusters 
(e.g., Richstone, Loeb, & Turner 1992). As we will see 
in §||, our model not only predicts fully developed ro- 
tation in A3266, but also predicts a viewing geometry 
that places a significant component of the opposing 
bulk flows along the observer's LOS, thus enhancing 
the prospect of detecting it spectroscopically. 

3. Observing the ICM Rotation with Astro-E 

The Astro-E XRS (Audley et al. 1999) is a high 
resolution X-ray spectrometer scheduled for launch 
in early 2000. With ~10 eV energy resolution and a 
quantum efficiency near unity across the 0.4 to 10 keV 
energy band, it represents the first opportunity to di- 
rectly observe gas dynamics in the ICM. Of course, 
these observations will be limited to LOS gas veloc- 
ities, and will be confused by multiple temperature, 
velocity, and possibly metal abundance components 
along the LOS. For this reason, it is extremely useful 
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to have a 3-dimensional model of the gas density, tem- 
perature, and dynamical structure both when plan- 
ning observations and interpreting the data. 

Figure || represents a simple LOS emission- weighted 
mean velocity map based on the numerical model. 
Analogous to the temperature map (Fig. |^), veloci- 
ties within a given computational zone along the LOS 
are simply weighted by the X-ray volume emissivity 
(oc n 2 T2). The range in velocities across the cluster is 
greater than 800 km s _1 . The opposing bulk motions 
evident on either side of the X-ray core constitute a 
strong signature of ICM rotation. West of the core, 
gas is moving away from the observer at greater than 
500 km s^ 1 while east of the core, gas is moving to- 
ward the observer at greater than 300 km s _1 . Of 
course, this velocity map is only indicative of the ob- 
servable LOS gas motions in the cluster. Below, we 
discuss detailed simulations of Astro-E observations. 

Based on the model, we propose two pointings se- 
lected to optimize their expected velocity separation 
and XSB. We suggest one pointing on each side of 
of the X-ray core, separated by ~8' (~ 600/iyo) on 
an east-west line running 2' north of the X-ray maxi- 
mum. The pointing location is chosen in part to avoid 
the point source located ~4' due east of the core. The 
expected mean emissivity- weighted velocities at these 
two locations are -195 km s _1 (east) and +533 km 
s" 1 (west) for a AV=728 km s _1 . 

In order to test the feasibility of this observation, 
we have used our numerical model and XSPEC (Ar- 
naud 1996) to generate synthetic X-ray spectra for 
each of the two pointings described above. Each 
spectrum is a composite of 50 spectra (Raymond- 
Smith+absorption) characterized by the local ICM 
temperature, density, velocity, and chemical abun- 
dance within a computational zone along the ob- 
server's LOS. We have assumed a uniform abun- 
dance of 0.2 solar (De Grandi & Molendi 1999). 
The total flux of the simulated cluster is scaled to 
2.84 xlO" 11 erg s" 1 cm" 2 (0.5-2.0 keV; David et 
al. 1999) while the absorption is characterized by 
N^.OxlO 20 cm" 2 (White et al. 1999). De Grandi 
& Molendi (1999) quote a somewhat lower absorp- 
tion of N h =1.6xl0 20 cm- 2 . Using the ROSATYSPC 
image as a flux distribution template, a list of pho- 
ton events is generated using MKPHLIST0. Optimum 
exposures are determined to be 50 ksec (east) and 60 
ksec (west). The resulting events were then input to 

■"-http: / /heasarc. gsfc.nasa.gov/docs/astroc/ 



XRSSIM 1 . After excluding photons not in the XRS 
FOV, the count rates for the east and west pointings 
are 0.20 and 0.17 cts s _1 , respectively, or approxi- 
mately 10 4 counts per spectrum. After rebinning of 
the photons to improve the statistics, we produce the 
spectra shown in Fig. |?]. An isothermal fit to the Fe 
K-line complex (6.2-6.7 keV) yields a velocity centroid 
shift of -219±85 km s _1 east of the core and +544±80 
km s _1 west of the core for a velocity separation of 
763±117 km s _1 between the two pointings. These 
values are consistent with those expected from a di- 
rect examination of the numerical model (see above) , 
thus demonstrating that the ICM rotation present in 
the model when scaled to the XSB of A3266 is in 
fact observable at a high level of significance with the 
Astro-E XRS. 

4. Summary 

We have presented a 3-dimensional numerical Hy- 
dro +N-body model of A3266 which is consistent 
with a wide range of observed properties. We believe 
A3266 represents an old (~3 Gyr), off-axis merger 
that is occurring into the plane of the sky at a 45° 
angle. The model is consistent, within resolution lim- 
its, with the ROSAT PSPC image, current ASCA 
temperature maps, BeppoSAX hard X-ray flux lim- 
its, the galaxy spatial and velocity distributions, and 
the existing radio data. We have also checked that the 
projected mass distribution agrees with gravitational 
lensing data. 

In this model, the off-axis merger has imparted sig- 
nificant angular momentum on the ICM of the merger 
remnant which we predict should be observable with 
the XRS on Astro-E. The signature of rotation will 
appear as two opposing bulk flows (AV~800km s _1 ) 
located on either side of the cluster's X-ray core along 
a line of constant declination. The degree of circula- 
tion present in the core of our model is consistent with 
high resolution Hydro+N-body simulations of galaxy 
cluster formation from cosmological initial conditions 
(Norman & Bryan 1998). Neither the galaxy redshift 
data nor the N-body particle distribution show simi- 
lar evidence of rotation. Unlike the ICM interaction 
which can be characterized as 'collisional', the inter- 
action between the subcluster galaxy components is 
'collisionless'. Therefore, while the transfer of angu- 
lar momentum between ICM components is very effi- 
cient, angular momentum is not transferred between 
the galaxy components. 
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It is difficult to assess the overall uniqueness of our 
model at the present time. Individually, no single 
observation places a strong constraint on the model. 
Even taken together, there is considerable flexibility 
in the merger parameters and viewing geometry, and 
we cannot rule out the possibility that recent mergers 
with significantly less massive subclusters have played 
some role (HDD99). Although limitations in the cur- 
rent data set certainly contribute to this problem, it 
is also indicative of an older merger that no distinc- 
tive features currently exist. In the event ICM rota- 
tion is not detected, the observations described here 
will nonetheless provide important new constraints on 
the model and on ICM dynamics in general. In ad- 
dition, these data will result in highly accurate tem- 
perature and metalicity measurements for two widely 
separated regions within a single cluster thus giving 
clues to the degree of spatial variation in both quan- 
tities. Inhomogeneities in the distribution of metals 
may help quantify the rate of mixing between gas 
components within the merger. 
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Fig. 1.— The ROSAT PSPC archival image of Abell 
3266. The image is background subtracted. Contour 
levels are 0.02, 0.04, 0.07, 0.14, 0.25, 0.45, 0.55, 0.65, 
0.75, 0.85, 0.95 of peak. 

Fig. 2. — A synthetic X-ray surface brightness image 
generated from the numerical model of A3266. Con- 
tours and linear dimensions (assuming H Q =70 km 
s _1 Mpc -1 ) are the same as in Fig. [l| The o's and 
+'s represent primary and secondary cluster particles, 
respectively. 

Fig. 3. — The projected emission-weighted temper- 
ature map of the A3266 model (color) overlaid with 
the X-ray surface brightness (contours). The num- 
bered boxes are the temperature regions defined by 
HDD99. 

Fig. 4. — Region-by-Region comparisons between the 
model temperature distribution and that observed by 
a) MFSV99 and b) HDD99 using ASCA. The regions 
defined by HDD99 are shown in Fig. [| 

Fig. 5. — Velocity vectors overlaid with X-ray emis- 
sivity (contours) within a 2-dimensional east-west 
(left to right) slice taken through the cluster core and 
along the observer's LOS (bottom to top). Note the 
large-scale counterclockwise rotation near the clus- 
ter core. Panel dimensions are 3.3 x 3.3 Mpc. The 
longest vector is ^900 km s _1 . 

Fig. 6. — Emission-weighted LOS ICM velocities 
(color) overlaid with X-ray surface brightness (con- 
tours). Velocities on the eastern half of the cluster 
(blue) are moving toward the observer while veloci- 
ties on the western half of the cluster (red) are moving 
away from the observer. The background color (light 
blue) corresponds to km s _1 . 

Fig. 7. — Simulated Astro-E spectra taken 4' on 
either side the X-ray core along an East- West line 
2' north of the core. Each represents a 50-component 
model (p,T,v from 50 zones within the numerical sim- 
ulation) normalized to the flux at the corresponding 
location in A3266. The vertical solid lines indicate 
the location of the fit to the line centroid. Integration 
times are 50 ksec East and 60 ksec West. 
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